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Abstract. We explore the connections between the well-studied Aztec Dia- 
mond graphs and a new family of graphs called the Half-Hexagons, discovered 
by Jonathan Novak. In particular, both families of graphs have very simple 
domino shuffling algorithms, which turn out to be intimately related. This con- 
nection allows us to prove an "arctic parabola" theorem for the Half-Hexagons 
as a corollary of the Arctic Circle theorem for the Aztec Diamond. 



In their groundbreaking paper |GV85j gave a method for counting families of 
non-intersecting lattice paths between two equinumerous sets of points. In their 
first example, the paths are between 



and are composed of unit-length steps up or to the right only. The number of such 
lattice paths is a Vandermonde determinant: 



Indeed, these ideas were already present in the works of Lindstrom |Lin73j and 
Karlin-McGregor [KM59j ; they are applicable far beyond the scope of these papers 
and form the enumerative-combinatorial cornerstone for many areas of modern 
mathematics. 

We shall focus on a specific case of the above example of Gessel-Viennot. Jonathan 
Novak pointed out to us that when Xi = 2i, then the number of these paths is 
2"( n + 1 )/ 2 ) which is the same as the number of domino tilings of an Aztec dia- 
mond [EKLP92aj; he asked us for a bijection. We didn't find one, but we did find 
many amazing similarities between these two models. Namely, they have similar 
domino shuffles, and similar limit laws. A little further detective work turned up 
a family of subgraphs of the Aztec diamond, which we call Aztec half- diamonds, 
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Figure 1. Five models which are in bijection: Non-intersecting 
lattice paths, pile of boxes, perfect matching and dual tiling, inter- 
lacing particle process, Staircase semistandard Young Tableau. All 
of these are order 3. 



Figure 2. Order 100 half- hexagon, as an interlacing particle process 

whose domino shuffling algorithm is identical, in a certain sense, to that on the 
half-hexagon. Our Aztec half-diamonds are similar but not identical to the half 
Aztec diamond of |FFllj . 

We wish to thank Jonathan Novak for bringing this problem to our attention, as 
well as our colleagues Alexei Borodin, Dan Romik, for helpful conversations. This 
paper began at the 2010 program in Random Matrix Theory, Integrable Systems 
and Interacting Particle Processes at the Mathematical Sciences Research Institute 
in Berkeley, California; it was completed while B. Young was visiting Universitat 
Wien. 

Several similar problems have been considered in the vast and ever-growing 
literature on the dimer problem, nonintersecting lattice paths, and the like. We 
mention some of them here. 

• Okounkov and Kenyon |KO07] calculate limiting shapes for dimer models 
on portions of the hexagonal grid, under "polygonal" boundary conditions; 
they show that for a generic polygonal boundary with 3d sides in which 
the edges appear in cyclic order, the limiting shape is an algebraic curve 
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of degree d. They comment that some of these conditions can be relaxed; 
however it is not clear to us how to use these methods to handle the erratic 
bottom boundary of our half-hexagons. 

• Di Francesco and Reshetikin [DFR09J study similar half-hexagonal shapes, 
but in which the long boundary is free; they obtain a variety of different 
limit shapes, none of which is the same as ours. 

• Borodin and Ferrari [BF08J have a very general framework for studying 
dynamics on interlacing particle processes, including the Aztec diamond 
domino shuffle and many others. Though they do not handle this particular 
case, our model does fit into this framework; we shall hopefully carry out 
their analysis in a subsequent paper. 

2. Bijective Combinatorics 

The Gessel-Viennot lattice path model is in bijection with a number of other 
combinatorial structures. Some of these bijections are "folklore" and all are well- 
known, but it is important to state briefly what they are, in order to establish 
terminology. 

For the remainder of this section, fix the order n e Z>o- 

2.1. Non-intersecting lattice paths. Let NILP(rz) be the set of families of non- 
intersecting lattice paths which begin at the points (0, — 2i) and end at the points 
(i, —i), composed of steps of unit length in the directions of increasing x and y. 
We get this from the example of [GV85J by taking X{ = 2%. The nth path from the 
top is of length 2n. 

As mentioned above, it is easy to enumerate these families of paths using the 
method of Gessel-Viennot; we shall do this in Section 13.11 



2.2. Lozenge tilings. It is well known (see, for example, |GV89l Poh05b] ) . that a 
family of non-intersecting lattice paths on the square lattice, with fixed start and 
end points, is in bijection with lozenge tilings of a certain region of the triangular 
lattice. Here, a lozenge is a parallelogram composed of two adjacent equilateral 
triangles. The boundary of the region depends only upon the locations of the 
endpoints. 

Definition 2.1. The regular triangular lattice C is the infinite planar graph whose 
vertices are the integer span of the vectors 

' 1 



1 

W= 2 



and which has edges joining any two vertices which are unit distance apart. L 
subdivides IR 2 into unit equilateral triangles. 

Let R n C IR 2 be the union of the large trapezoid with corners 

{nv, nw, n(w — v), —nv} 
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0, 



n — 1 



U7, 



w. 



and the n small trapezoids with corners 

(2i — n + l)v + v, 
(2i — n + l)v — v, 
(2i — n + l)v — w, 
(2i — n + l)v + v — w 

See Figure [IJ pictures 2 and 3, for a graphical representation of R3. Let LT(n) 
be the set of lozenge tilings of R n . To obtain an element of LT(n) given an element 
of NILP(n), first apply the affine transformation which takes the points 

(0, -2i) 1 — > {-n + 2i)v- \ 

{i, —i) 1 — > -nv + (i — |) 

Now the lattice steps are in directions — v and w — v; each step begins and ends 
on the boundary of a triangle in C and traverses two of the triangles of L. Form 
a partial tiling by placing the corresponding lozenge over each step. The holes in 
this tiling can be covered uniquely by lozenges as well. 

Observe that, in drawing this tiling, we have also drawn a pile of cubical boxes: 
each tile represents a visible face of a cube with faces parallel to the coordinate 
planes, viewed isometrically from the direction (1,1,1). 

2.3. Perfect matchings on a half-hexagon. Let £ v denote the planar dual 
of C. £ v is the regular tiling of the plane with hexagons, sometimes called the 
honeycomb mesh, grid or lattice. 

Definition 2.2. Let be the subgraph of C v induced by those vertices which lie 
within R. We call R^ the half-hexagon graph. Let HH{n) denote the set of perfect 
matchings of R^ . 

There is a folklore bijection between HH(n) and LT(ra). Suppose we are given 
an element T of LT(n). Each lozenge in the tiling T is composed of two triangles 
of C, which are dual to two vertices in £ v . Join every such pair of vertices with an 
edge to obtain a matching in HH(ra). The resulting matching is perfect (i.e. that 
every vertex is covered) by virtue of the fact that a tiling in LT(n) covers R(n) 
completely. 



2.4. Interlacing particle process. Observe that the vertical edges of 
centered at the points 



arc 



-n — - + j ] v + [n — i)w 



n; j = 1, . . . , i + n 



In the above set, we call i the row index and j the position. Let tt g HH(n). 
Observe that the vertical edges in tt determine n completely, subject to the fol- 
lowing interlacing condition: if there are edges in positions j and j' of row i, then 
there must be an edge in some position j" of row i — 1, with j < j" < j. Indeed, 
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any collection of vertical edges which interlace in this manner determine a tiling 
in HH(rt). 

From 7i~, then, we may construct an interlacing particle process 

| (—n — \ + j)v + (n — i)w is the center of a vertical edge in n } C N 2 . 
Let IPP(n) denote the set of all such interlacing particle processes. 

2.5. Staircase tableaux. Let P be an element of I PP(n). Define the numbers gij 
by 

n-l 

P={J{(hg lJ )\o<3 <t-i} 

i=0 

In other words, is the position of the jth vertical edge in row i of the corre- 
sponding perfect matching in HH(n). The interlacing conditions imply that 

9ij < 9i-l,j < 9i,j+li 

in other words, that (g^) is a Semistandard Young tableau of staircase shape |Sta99j 
with bottom row equal to (1, 3, 5, . . . , 2n + 1). We call these objects Staircase 
tableaux for short, and we write the set of all such as ST(n). Note also that the 
numbers 

form a Gelfand-Tsetlin pattern (see |Sta99| (7.37)]) with bottom row equal to 
(0,1,2,. ..,n). 

3. Enumerative Combinatorics 

There are at least two easy ways to enumerate HH(n) by evaluating determi- 
nants. We shall introduce a third way as a consequence of the shuffling algorithm. 
Doubtless there are many others. 

3.1. Non-intersecting lattice path enumeration. Gessel-Viennot [GV85] han- 
dle a slightly more general situation. Observe that there are ("?) up-right lat- 
tice paths from the point (0, — Xj) to (j, — j). Consequently, given n integers 
X\ < X2 < ■ ■ ■ < x n , the number of families of nonintersecting lattice paths 
from {(0, — Xi] to {(j, — j)} is given by 0. This determinant evaluation is given 
in |GV85] : it works because it is essentially a Vandermonde determinant. Putting 
Xi = 2i, we recover the endpoints of the paths for NILP(n), and we obtain 



NILP(n)| = 2 n(n+1)/2 . 
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3.2. Staircase tableau enumeration. An alternate easy enumeration, this time 
of ST(n), goes through symmetric function theory. First we observe |Sta99l (7.37)] 
that because of the Gelfand-Tsetlin pattern interpretation, we have 

|ST(n)| = s A (l,...,l) 

where s\ is a Schur function with n+1 arguments and A is the "staircase partition" 
(1, 2, • ■ ■ , n + 1). It is a consequence of the classical bialternant definition of the 
Schur function (see [Sta991 Chapter 7.15 and ex. 7.30] that 

S X {xi,. . . ,X„+i) = [ J (Xi + Xj). 

\<i<j<n 

Putting all Xi equal to 1 gives |ST(n)| = 2 n( - n+1 ^ 2 as before. Indeed, performing 
the principal specialization Xi i— y q l gives a g-enumeration of HH(n), in which each 
element tt of HH(n) is assigned weight proportional to q vol ( n \ Here, vol(7r) denotes 
the integral of the height function of tt. See |Ken04j for an introduction to the 
height functions of dimer models. 

4. Dynamics 

4.1. The half-hexagon shuffle. In this section, we will define dynamics on in- 
terlacing particles, called the half-hexagon domino shuffle, which takes the form 
of a random map from HH(n) to HH(n + 1). The procedure is easiest to describe 
and implement on the staircase tableaux ST(n), and it is easiest to illustrate on 
the half-hexagon perfect matchings HH(n) (see Figure [3]). 

The half-hexagon domino shuffle constructs a perfect matching tt' on R v (n + 1) 
from tt. It is not a deterministic process: one has to make a series of fair coin tosses 
to determine tt' . As it turns out, these coin tosses provide a direct explanation of 
the fact that the cardinality of HH(n) is a power of two. 

Algorithm 4.1. (Domino shuffling for the half-hexagon) 
Input: 

(dij), o staircase tableau in ST(n), 

independent Bernoulli 0-1 variables. 

Output: 

(hij), a staircase tableau in ST(n + 1). 

1: for i from to n: 

f or j from to i: 

3: if j < i and g^j = hi-xj then: 

4'- hij gij 

5: else if j > and gij = h^xj-i then: 

6: h id <r- g itj + 1 

7: else 
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Figure 3. Domino shuffling on the half-hexagon. Groups of edges 
which are forced to move together, in order to maintain the inter- 
lacing conditions, are circled. 



8: hij «— g i:j + £y 

9'- for j from to n: 
10: h n+1J 2j • 1 

It is easier to illustrate this algorithm acting on the half-hexagon (see Figure [3]) 
though slightly harder to describe. Begin with a perfect matching it on HH(n) 
(specified by its vertical edges). Observe that it is possible to add, deterministi- 
cally, a row of n + 1 vertical edges, appearing every second edge, to the bottom of 
£ v , in such a way as to maintain the interlacing condition. Further, it is possible 
to superimpose the next-larger half-hexagon i? v (n + 1) on this graph, in such a 
way that each vertical edge is in the center of a hexagon of R v (n + 1). 

Working from the top of the graph to the bottom, each vertical edge in ir jumps 
either left or right onto the nearest vertical edge in the same row of R v (n + 1 ) . 
These jumps happen independently at random, according to the result of a fair 
coin toss, with the following exceptions: 

• if moving an edge left would violate the interlacing condition with the new 
row above, then the edge moves right with probability 1. 

• if moving an edge right would violate the interlacing condition with the 
new row above, then the edge moves left with probability 1. 

We invite the reader to check the preceding procedure is the same as Algo- 
rithm 14.11 



One does need to check that the output of Algorithm |4.1| is always a staircase 
tableau, i.e. to verify that the interlacing conditions hold. This is done inductively 
on the row i, together with a checking that the deterministic row n + 1 interlaces 
with row n. 
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4.2. Preserving the uniform distribution. In this section, we argue that the 
domino shuffle described above preserves the uniform distribution (using the ter- 



minology of building rules). We first introduce the time reversal of Algorithm 4.1 



Algorithm 4.2. (Time-reversed domino shuffle) (Time-reversed domino shuffling 
for the half-hexagon) 

Input: 

(gij), a staircase tableau in ST(n + 1), 
(iij), independent Bernoulli 0-1 variables. 

Output: 

(hij), a staircase tableau in ST(n). 



11 

12: 
13: 

H 
15: 
16: 
17: 
18: 
19: 
20: 



for j from to n: 

h nd = 2j + 1 
for i from n — 1 down to : 

f or j from to i: 

ifgij = h i+1J then: 

^ Qi,j 

else if gij = /ii+ij+i then: 

hi,j ^ Qi,j 1 

else 

^ 9i,j 



On the half hexagon, this algorithm does the following: the (n + l)st row of 
vertical edges is dropped altogether; the nth row jumps deterministically to posi- 
tions 0, 2, . . . , n. Then, working from the bottom to the top, edges jump left or 
right with probability |, except that edges are sometimes forced to move in order 
to interlace with edges in the row below. 

Definition 4.3. Let CHH(n) denote the vector space whose orthonormal basis is 
indexed by the elements of HH(n). Let (•, -) n denote the inner product which makes 
this an orthonormal basis. 

Let P(tt — > 7r') be the probability that Algorithm 4-1 produces output tt' when 
given input n. 

Let ip '■ CHH(n) — > CHH(n + 1) be the linear map for which 

(V>(7r),7r') = P(vr ->- tt'). 

Lemma 4.4. The maps ip an d i>' are adjoint to each other. That is, if tt £ 
HH(n - 1) and ip G HH(n), then 
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Proof. Algorithms 4.1 and 4.2 move the edges either randomly (steps [8] and 20) or 
deterministically (steps El [6 16 and 18). We call the random moves free and the 



deterministic moves forced, because all of the deterministic moves occur in blocks, 
precipitated by a preceding free move. Figure [3] shows an instance of the shuffling 
algorithm, with blocks of forced moves circled. In this terminology, we have 

ib 

^ n _ 2-#{free choices in tt -> tt'} 

ib' 

1^ ^'tjA — 2~#{ free choices in tt' -> tt} 

Let S be all of the free choices in tt — > tt ; and let S be the set of all free choices 
which do not force any edges in row n of tt to move. Since there are n edges in 
the bottom row, \S\ = \S'\ + n. 

Let ej be one of the edges in E in row i, and suppose that in passing from tt 
to tt'. Say that forces the movement of edges e i+ i, . . . , e#, in successive rows. 
Observe, then, that the time-reversed algorithm if}', in passing from tt' to tt, sees a 
free choice at edge ey which forces edges e«, . . . , ei + \ to move. Moreover, tt' moves 
row n deterministically (and then deletes it) so these are never free choices. As 

such, the free choices in tt — > tt are in bijection with 5". 

(iPtt,tt') = 2~ |s| = 2' 15 ' 1 "" = 2- n (vr,^V) 

□ 

Definition 4.5. Let fi n G CHH(n) denote the vector corresponding to the uniform 
probability distribution on HH(n): 

^ n= \HH(n)\ n - 

Proposition 4.6. Domino shuffling preserves the uniform distribution: 

tf>Hn-i = b-n- 
Proof. This is a straightforward computation: 

TreHH(n-l) Tr'eHH(n) 



2 Ti 

Tr'eHH(n) \vreHH(n-l) 

= - E - 

2 n J 

Tr'£HH(n) 

In the second line, the bracketed sum is equal to one for any tt' because if>' is 
stochastic. The latter vector is proportional to /i n , and is thus equal to \x n because 
the coordinates of both vectors sum to one. □ 
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Note that it follows from this proof that \HH{n)\ = 2 n \HH(n — 1)|, which 
provides a new derivation of the fact that \HH(n)\ = 2™( n+1 )/ 2 . 

5. Review: the Aztec diamond 

The Aztec Diamond of order n (which shall be denoted A n ) is a certain shape 
in the plane that can be covered with dominoes in 2™( n+1 )/ 2 ways. More precisely, 
A n consists of the union of all squares whose corners have integer coordinates, 
whose sides are parallel to the coordinate axes and whose interiors are contained 
in {(x,y) G IR 2 : \x\ + \y\ < n + 1}. The model was introduced in |EKLP 92a, 
IEKLP92bj in the study of Alternating sign matrices and a good survey is |Joh05a] . 
In this section we give a short review of the results we need about this model with 
citations. 

One way to study the random tilings in this model is to introduce a certain 
particle process, see |Joh05a[ INorlO] and also Figure [5] for details. In short, the 
lattice on which the dominos are placed is colored like a chessboard, and all domi- 
noes whose bottommost or rightmost square is dark is represented by a particle. 
The tiling does uniquely determine the positions of the particles. The converse is 
only true modulo the fact that there will be 2 x 2-rectangles where two horizontal 
dominoes can be switched for two vertical ones or the other way around. Particle 
configurations coming from a tiling of A n are in bijection with Alternating sign 
matrices and this is the original motivation for studying this tiling model. How- 
ever, the measure on particle configurations that is induced by uniform measure 
on all possible tilings is not the same as uniform measure on all Alternating sign 
matrices (it is, rather, connected to the 2-enumeration of alternating sign matrices, 
see |EKLP92aj ). 

In a tiling of A n there are particles on n rows, see Figure |4j Along the 

line y — 1 there is a single particle, along line 2 there are two, etc. Let Xj be the 
position of the jth particle on line y — i. Observe that the particles interlace in 
the sense that 

xl +1 <x\< x{~l] 

this is similar, but not the same, as the interlacing conditions for the half-hexagon 



particle process in Section 2.4 



There is an algorithm, called the shuffling algorithm, that can be used to con- 
struct a random tiling and is described in great depth in each of |EKLP92bjl IJPS981 
IPro03] . The procedure starts with a tiling of A n ; one moves the dominos of the 
tiling about in a certain way, decided by a certain number of coin flips, producing 
a tiling of A n+ i. If one starts with the empty tiling of A Q , performs this algorithm 
n times, tossing fair coins all the way, then one ends up with a sample from the 
uniform distribution of all possible tilings. 

The idea of |NorlO] is to look at the positions of the aforementioned particles 
under the evolution of this algorithm. It turns out that the particles are not glued 
to the tiles. The dynamics of these particles are as follows. In the following all 



DOMINO SHUFFLING ON NOVAK HALF-HEXAGONS AND AZTEC HALF-DIAMONDS 11 




Figure 4. Particle processes associated with domino tilings of an 
Aztec diamond. The underlying lattice is colored like a chessboard 
and all dominoes whose bottommost or rightmost square is dark is 
represented by a particle. 

7*(i) for i, j, t = 1, 2, . . . , are independent Bernoulli random variables, that is 
one with probability \ and zero otherwise. 

It turns out that the first particle performs the simple random walk 

(2) x\(t)=x 1 1 (t-l)+ 1 1 1 (t). 

The particle x\ performs a simple random walk with a reflecting boundary. More 
precisely, while x\(t) < x\(t) it performs a random walk independently of x\, at 
each time either staying or adding one with equal probability. However, when 
there is equality, x\(t) = x\(t), it is pushed forward by that particle. In order to 
represent this as a formula, we subtract one if the particle attempts to jump past 
x\. 

(3) xl(t) = x\{t - 1) + 7l 2 (t) - l{x{(t - 1) + 7l 2 (i) = x\(t-l) + l} 
By symmetry then 

(4) x\{t) = x\(t - 1) + 72 2 (t) + \{x\(t - 1) + 72 2 (t) = x\{t- 1)}. 
The same pattern repeats itself evermore. 

(5) x{{t) = 4(t - 1) + 7 ?'(t) - l{x{(t - 1) + 7 ?'(t) = x{~\t - 1) + 1} 

(6) x]{t) = ,••;(/ - 1) + 7j(t) + - 1) + 7j(t) = :(/ - 1)} 

(7) x>{t) = x{(t - 1) + T>(t) - l{4(t - 1) + l({t) = z>-\t - 1) + 1} 

(8) +l{xi(t-l)+ 1 {(t)=xp L (t-l)}. 

with initial conditions = i for 1 < i < j. 

In the analysis [NorlOj of the asymptotics of the domino shuffling algorithm, it 
turns out to be quite inconvenient that the n particles on level n are not created 
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Figure 5. The height function is defined on the lattice on which 
the corners of the dominoes line up. 



until time n. A simple change of variables will fix this. Let 

(9) Xl(t)=x\(t + i) 

for 1 < % < j and t — 1, 2, ... . We mention this here for a completely different 
reason: Rewriting the above recursion formulas in terms of the variables (A/)!<j<,,- 



gives the same recursion as is implemented by Algorithm 4.1 above (though with 



a different initial condition, see Section 6.1). 



6. The Aztec Half-Diamond 

Recall the definition of the Aztec diamond A n in Section [5j The Aztec half- 
diamond H n of order n is a certain subregion of A n . More precisely, for n even, 

H n = {(x,y)EA n :y>2[x/2\} 

and, for n odd, 

H n = {{x,y) G A n : y + 1 > 2[(x + 1)/2J}. 

Though the definition as stated is a bit cryptic, the Figure [6] should make this 
quite clear. 

6.1. Dynamics: Domino shuffling. The Aztec Half-Diamonds also have a domino 
shuffling algorithm. Propp |Pro03] describes a way to generate random perfect 
matchings on certain planar bipartite graphs called "generalized domino shuffling" . 
The strategy is to embed the graph 67 into an Aztec diamond of sufficiently large 
order n, and then compute a certain series of probabilities Pj(m) where 1 < m < n 
and 1 < i,j < m (this is perhaps an overly concise summary, but both the man- 
ner of embedding and the means of computation are described explicitly enough 
in |Pro03j to allow computer implementation). To generate the perfect matching, 
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Figure 6. Two Aztec half-diamonds make an Aztec diamond. The 
left picture shows half-diamonds of orders 18 and 20; the right pic- 
ture orders 19 and 21. The interlacing particle process is also shown. 
The right picture is the domino shuffle of the left. 



one then simply runs the domino shuffling algorithm, with the following modi- 
fication: make Q(m) a Bernoulli random variable which takes the value 1 with 
probability P-(m). The final tiling of the Aztec Diamond of order n, restricted 
to the embedded copy of G, turns out to be uniformly random; the intermediary 
tilings of the smaller Aztec diamonds may be discarded. 

Using the embedding shown in Figure [6] for generalized domino shuffling, one 
can check that the intermediary tilings of the smaller Aztec diamonds are also of 
the form of Figure [6] As such, Propp's algorithm is a "domino shuffle" for the 
Aztec half-diamond in our sense: a random, locally defined map which increases 
the size of the half-diamond but preserves the uniform distribution. Indeed, this 
procedure is exactly the same as ordinary domino shuffling everywhere except the 
center line when n is even; at these times, the particles in the center are forced to 
jump to equally spaced positions. This is equivalent to imposing condition 

x™ +1 (2m + 1) = 2i 

for 1 < i < m. 

6.2. Height functions and limit shape. It is possible |Thu90] to associate a 
discrete surface in IR 3 , called a height function, to any domino tiling of the plane. 
More precisely, the height function h : Z 2 — > Z, where the domain is the square 
grid whose vertices coincide with the corners of the dominos and the centers of 
their edges. 
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Definition 6.1. Let T be a domino tiling of the region R. Then h : Z 2 — > Z is a 
height function for T if, whenever x + y = (mod 2) and (x, y) is in R, then 

• h(x, y) = h(x, y + 1) + 1 if the edge from (x, y) to (x, y + 1) does not cross 
a domino in T; 

• h(x, y) = h(x + 1, y) — 1 t/ie edge from (x, y) to (x + 1, y) does not cross 
a domino in T. 



Note that h determines T uniquely, and that two height functions for T differ 
only by a constant. 

This definition coincides with those in |Thu90l ICKPOlj and, In the case where 
the region R is an Aztec Diamond, this definition appears in [E KLP92aj . where it 
is closely related to the height function for an alternating sign matrix. 

The reader should be advised that there are closely related concepts in the liter- 
ature called relative height functions |KO07t IKOS06] , edge-placement probabilities 
and one-point functions for the particle process |Ken97t [Pro03| . 

Fix a region R, a tiling T of R, and a height function h for T. Since no tiles ever 
cross the boundary of R, the restriction of h to OR is independent of T. Indeed, 
h\ m can even be computed without specifying T at all. 



Definition 6.2. The function h \qr is called a boundary height function. 



Proposition 6.3. (See |Thu901 Section 4]j R possesses a domino tiling if and 
only if R has a well-defined boundary height function. 



Once it became possible to generate uniformly random tilings of large Aztec dia- 
monds, it became immediately obvious that all such tilings have the same "shape" . 
It was first shown in |JPS98] that a typical tiling of a large Aztec diamond has all 
of its disorder concentrated in a circular region, where the circle is tangent to all 
four sides of the diamond; the tiling is frozen close to the four corners. 

In fact, more is true: the height function of a random tiling of an Aztec diamond 
tends to the following explicit limit. This limit seems possible to do using the 
correlation kernel for the Aztec Diamond, as defined in[Joh05aJ, and something 
similar (asymptotics for edge placement probabilities) were computed in [CEP96J), 
but the first explicit derivation would seem to be in |Rom09j . The coordinates 
are slightly different: in the following theorem, h*- represents a rescaled height 
function of an order n Aztec diamond, where the domain is to [0, 1] x [0, 1] and the 
range is rescaled to [0, 1]. Indeed, for the remainder of this section we will work in 
these coordinates. 
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Theorem 6.4. (Theorem IV in |Rom09] ) Define 
Z(x,y) 



2 



(x — 1/2) arctan 




I {x _ l/ 2 )2 _ {y _ l/2) 2 



1/2 -y 



+ - arctan 
2 



2Qr-l/2)(l/2-y) 

v/| - - I/ 2 ) 2 - (y - 1/2) 2 



— (1/2 — w) arctan 




x 



1/2 



| - (x - 1/2) 2 - (2/ - 1/2)= 



and define 



G[x,y] 



x + y 

x + Z(x, y) 
x-y 



0<z<^4^, 

1-20,(1-;/) 1+20(1-1/) 



l+20(l-y) 



< X < 1. 



Tnen as n — )■ oo we nave iae convergence in probability 

G(i/n,j/n 



max 

0<i,j<n 



h* 1 



n 



> 0. 



(10) 



G[x, 1/2] 



In comments after Equation (15), |Rom09j observes that 

1 

2 ! 

that is to say, that the limiting height function is constant across the center of the 
Aztec Diamond. (It is also possible to observe this directly: the Aztec diamond 
has a reflection symmetry in the line y = 1/2 which leaves the uniform measure 
invariant, but which negates height functions up to a constant). This observation is 
quite important for our purposes: it means that the limit shape for Aztec Diamonds 
coincides with the boundary height function for Aztec Half-Diamonds. 

If we replace the Aztec diamonds with a sequence of regions R n approximating 
a given region R, in such a way that the rescaled boundary height functions also 
tend to a limiting boundary height function on R, then there is still a unique 
limiting shape and a variational principle (though in general there is no reason to 
expect this limiting shape to be as nice as a circle!) 

Theorem 6.5. ( |CKP01] Theorem 1.1] Let R* be a region in IR 2 bounded by a 
piecewise smooth, simple closed curve dR* . Let : OR* — > R be a function which 
can be extended to a function on R* with Lipschitz constant at most 2 in the sup 
norm. Let f : R* — > R be the unique such Lipschitz function maximizing the 
entropy functional Ent(f), subject to f\dR* = hb- 
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Let R be a lattice region that approximates R* when rescaled by a factor of 
1/n, and whose normalized boundary height function approximates h^. Then the 
normalized height function of a random tiling of R approximates f , with probability 
tending to 1 as n — >■ oo. 

In the above theorem, Ent(/) is given by 

and ent fj^J is a certain explicit function of the gradient of the surface h; 

see |CKP01] for further details. For us, the most important point is that the 
asymptotic number of tilings depends only upon the asymptotic height function 
g, not upon the precise nature of the boundary of R; moreover this dependence is 
local. As such, an immediate consequence of Theorem 6.5 is the following: 

Lemma 6.6. Let S C R be regions. Let hg R be a boundary height function on 
OR, and let f : R — » K be the entropy-maximizing asymptotic height function for 
(R, hg R ). Let g be the entropy-maximizing asymptotic height function for (S, f\as)- 
Then g = f\ s . 

Proof. We will write Ents(.), Entjj(.), Entjj\g(.) for the entropy functionals of the 
various regions, and As, A R , A r \ s for their areas. It is a direct consequence of (11 ) 
that 

A R Ent R (f) = AsEnt s (f) + A R \ s Ent R \ s (f)- 
Define a height function g* on R as follows: 



(12) g*(x) 




if x G S 
\txeR\S 



Observe that this function is a well-defined asymptotic height function because 
g = f on the boundary of S. If Eiats(f\ds)) < Ents(g), then Eiat R f < Eiat R g*, 
contradicting the fact that / is the entropy maximizer on R. Thus Ents{f\ds)) — 
Ent s{g), and so g = f\g is the unique entropy maximizer on S. □ 

Corollary 6.7 . Th e limit shape of the Aztec half-diamond is the restriction of G 
(from Theorem 6.4) to y < |. 

Proof. Let 

R= [0,1] x [0,1] f(x,y) = G[x,y] 

S=[o,i] x [0,1/2] g\ds = f\ds; 

That is, R and / describe the limit height function for the Aztec Diamond, and S 
is the limiting region for the Aztec Half-Diamond, and g\gs is its limiting height 
function by Equation (10). Lemma 6.6 now says that the limiting height function 
for the Aztec Half-Diamond g coincides with fs- □ 
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Figure 7. Space-time diagram for domino shuffling. Left image: 
the Aztec diamonds and their arctic circles. Right image: the half- 
hexagon and its arctic parabola. 

Corollary 6.8. The limit shape of the half-hexagon is a parabola. 

Proof. One checks that the particle process defined by the half-hexagon domino 
shuffle has the Aztec Half-Diamonds as its constant-time slices (see Figure [7]). 
Indeed, the behaviour of Algorithm 4.1 coincides with the equations |9j except 



that the initial condition is different; it straightforward to check, that this initial 



condition is enforced both by Algorithm 4.1 and by Propp's generalized domino 



shuffling [Pro03] applied to the half-hexagon; see Section 6.1 

The space-time diagram of domino shuffling traces out a cone inscribed in a 
square pyramid. The half-hexagons' perfect matchings correspond to slices parallel 
to one of the sides of this cone; the intersection is a parabola. 

□ 

We remark that one can in fact write down the limiting height function in this 



way as well; it is precisely the image of the height function of Theorem 6.4 under 
the affine transformation which takes the rectangle [0, 1] x [0, 1/2] to the trapezoid 
with corners {(±1,0), (±(1/2), y / 3/2)}. Moreover, domino shuffling on the Aztec 
half-diamond has a reasonably nice description: it turns out to coincide with a 
particular instance of Propp's Generalized Domino shuffling [Pro03j up to gauge 
transformation. 



7. Future work and open problems 

The reader will doubtless have noticed that this paper is a study of the phe- 
nomenology of one of the nicest, most special cases of the dimer model on planar 
bipartite graphs; naturally one should try to push the analysis to more general 
situations. For instance, a few other domino shuffles have been discovered on a 
variety of statistical mechanical models (see, for example, [YClOt IBG09] ) and it 
would be very interesting to study different space-time sections of them. 



18 



ERIC NORDENSTAM AND BENJAMIN YOUNG 



Taking slices through the cone in Figure [7] with varying slopes should yield a 
family of particle processes whose limit shapes are conic sections. Though this is a 
fairly trivial observation, we note that relatively few instances of the dimer model 
have known low-degree algebraic limit shapes. 

There is a general framework [BF08J for studying dynamics of the sort we have 
described in this paper, geared in particular to asymptotics. It appears that our 
model fits into this framework, but we have yet to determine whether the compu- 
tations are tractable. 

We conclude with the following remark: the problem which Jonathan Novak 
posed to us — give a bijection between AD(n) and HH(n) — is still open, and we 
would be very happy to see a solution! 
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